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Abstract 

In this paper we propose the time & fitness-dependent Hamiltonian form of human biome- 
chanics, in which total mechanical + biochemical energy is not conserved. Starting with the 
Covariant Force Law, we first develop autonomous Hamiltonian biomechanics. Then we extend 
it using a powerful geometrical machinery consisting of fibre bundles, jet manifolds, polysym- 
plectic geometry and Hamiltonian connections. In this way we derive time-dependent dissipa- 
tive Hamiltonian equations and the fitness evolution equation for the general time & fitness- 
dependent human biomechanical system. 

Keywords: Human biomechanics, configuration bundle, Hamiltonian connections, jet mani- 
folds, time & fitness-dependent dynamics 
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1 Introduction 



Most of dynamics in both classical and quantum physics is based on assumption of a total energy 
conservation (see, e.g. [4]). Dynamics based on this assumption of time- independent energy 
usually given by Hamiltonian (or Lagrangian) energy function, is called autonomous. This basic 
assumption is naturally inherited in human biomechanics, formally developed using Newton-Euler, 
Lagrangian or Hamiltonian formalisms (see [TJ ® 131 M M HI GO HO])- 

And this works fine for most individual movement simulations and predictions, in which the 
total human energy dissipations are insignificant. However, if we analyze a lOOm-dash sprinting 
motion, which is in case of top athletes finished under 10 s, we can recognize a significant slow- 
down after about 70 m in all athletes - despite of their strong intention to finish and win the 
race, which is an obvious sign of the total energy dissipation. This can be seen, for example, in 
a current record-braking speeddistance curve of Usain Bolt, the world-record holder with 9.69 s, 
or in a former record-braking speeddistance curve of Carl Lewis, the former world-record holder 
(and 9 time Olympic gold medalist) with 9.86 s (see Figure 3.7 in 10]). In other words, the total 
mechanical + biochemical energy of a sprinter cannot be conserved even for 10 s. So, if we want 
to develop a realistic model of intensive human motion that is longer than 7-8 s (not to speak for 
instance of a 4 hour tennis match), we necessarily need to use the more advanced formalism of 
time-dependent mechanics. 

In this paper, we will first develop the autonomous Hamiltonian biomechanics as a Hamiltonian 
representation of the covariant force law [see |T]) in the next section] and the corresponding covariant 
force functor. After that we will extend the autonomous Hamiltonian biomechanics into the time- 
dependent one, in which total mechanical + biochemical energy is not conserved, for this, we will 
use the modern geometric formalism of jet manifolds and bundles. 

2 The Covariant Force Law 

Autonomous Hamiltonian biomechanics (as well as autonomous Lagrangian biomechanics), based 
on the postulate of conservation of the total mechanical energy, can be derived from the covariant 
force law [H [2j El H] , which in 'plain English' states: 

Force 1-form = Mass distribution x Acceleration vector-field, 

and formally reads (using Einstein's summation convention over repeated indices): 

F t = m lj a j . (1) 

Here, the force 1-form Fi — Fi(t, q 7 p) = F!(t, q,q), (i = 1, n) denotes any type of torques and 
forces acting on a human skeleton, including excitation and contraction dynamics of muscular- 
actuators [13l [12l [11] and rotational dynamics of hybrid robot actuators, as well as (nonlinear) 
dissipative joint torques and forces and external stochastic perturbation torques and forces [5]. 
rriij is the material (mass-inertia) metric tensor, which gives the total mass distribution of the 
human body, by including all segmental masses and their individual inertia tensors. a J is the total 
acceleration vector-field, including all segmental vector-fields, defined as the absolute (Bianchi) 
derivative v l of all the segmental angular and linear velocities v l = i 1 , (i = 1, n), where n is 
the total number of active degrees of freedom (DOF) with local coordinates {x 1 ). 
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More formally, this central Law of biomechanics represents the covariant force functor 
defined by the commutative diagram: 



TTM 




(2) 



Here, M = M n = {x % , (i — l,...,n)} is the biomechanical configuration n— manifold, that 
is the set of all active DOF of the biomechanical system under consideration (in general, human 
skeleton), with local coordinates (x l ). 

The right-hand branch of the fundamental covariant force functor T„ : TT*M — ► TTM de- 
picted in @ is Lagrangian dynamics with its Riemannian geometry. To each n— dimensional (riD) 
smooth manifold M there is associated its 2nD velocity phase-space manifold, denoted by TM 
and called the tangent bundle of M. The original configuration manifold M is called the base of 
TM. There is an onto map tt : TM — ► M, called the projection. Above each point x G M there 
is a tangent space T X M = w~ 1 (x) to M at x, which is called a fibre. The fibre T X M C TM is 
the subset of TM, such that the total tangent bundle, TM = | | T X M, is a disjoint union of 

m£AI 

tangent spaces T X M to M for all points x G M . From dynamical perspective, the most important 
quantity in the tangent bundle concept is the smooth map v : M — * TM, which is an inverse to 
the projection n, i.e, n o v = Mm, tt(v(x)) = x. It is called the velocity vector-field v l = x l Its 
graph (x,v(x)) represents the cross-section of the tangent bundle TM. Velocity vector-fields are 
cross-sections of the tangent bundle. Biomechanical Lagrangian (that is, kinetic minus potential 
energy) is a natural energy function on the tangent bundle TM. The tangent bundle is itself a 
smooth manifold. It has its own tangent bundle, TTM. Cross-sections of the second tangent 
bundle TTM are the acceleration vector-fields. 

The left-hand branch of the fundamental covariant force functor T* : TT*M — > TTM depicted 
in ^ is Hamiltonian dynamics with its symplectic geometry. It takes place in the cotangent bundle 
T*M ro b, defined as follows. A dual notion to the tangent space T X M to a smooth manifold M at a 
point x = (x l ) with local is its cotangent space T*M at the same point x. Similarly to the tangent 
bundle TM, for any smooth nD manifold M, there is associated its 2nD momentum phase-space 
manifold, denoted by T*M and called the cotangent bundle. T*M is the disjoint union of all its 
cotangent spaces T*M at all points x G M, i.e., T*M = [J T*M. Therefore, the cotangent 

bundle of an n— manifold M is the vector bundle T*M — (TM)*, the (real) dual of the tangent 
bundle TM. Momentum 1 -forms (or, covector-fields) pi are cross-sections of the cotangent bundle. 
Biomechanical Hamiltonian (that is, kinetic plus potential energy) is a natural energy function on 

1 This explains the dynamical term velocity phase-space, given to the tangent bundle TM of the manifold M. 
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the cotangent bundle. The cotangent bundle T*M is itself a smooth manifold. It has its own 
tangent bundle, TT*M . Cross-sections of the mixed-second bundle TT*M are the force 1-forms 
Fi = Pi. 

There is a unique smooth map from the right-hand branch to the left-hand branch of the 
diagram 

TM 3 {x\v l ) h-> (x^p*) g T*M. 

It is called the Legendre transformation, or fiber derivative (for details see, e.g. [21 3]). 

The fundamental covariant force functor : TT*M — ► TTM states that the force 1-form 
Fi = pi, defined on the mixed tangent-cotangent bundle TT*M, causes the acceleration vector- 
field a 1 — v l , defined on the second tangent bundle TTM of the configuration manifold M. The 
corresponding contravariant acceleration functor is defined as its inverse map, T* : TTM — ► 
TT*M. 

Representation of human motion is rigorously defined in terms of Euclidean S'i?(3)-groupa3 of 
full rigid-body motion in all main human joints ,6[. The configuration manifold M for human 
musculoskeletal dynamics is defined as a Cartesian product of all included constrained SE(3) 
groups, M — Y[j SE(3y where j labels the active joints. The configuration manifold M is coordi- 
nated by local joint coordinates x l (t), i = 1, n = total number of active DOF. The corresponding 
joint velocities x l {t) live in the velocity phase space TM, which is the tangent bundle of the con- 
figuration manifold M. 

The velocity phase-space TM has the Riemannian geometry with the local metric form: 



(g) = ds 2 = gijdx l dx\ 



(3) 



where gij(x) is the material metric tensor defined by the biomechanical system's mass-inertia 
matrix and dx % are differentials of the local joint coordinates x % on M. Besides giving the local 
distances between the points on the manifold M, the Riemannian metric form (g) defines the 
system's kinetic energy: 

T = —gtjX x , 

giving the Lagrangian equations of the conservative skeleton motion with kinetic-minus-potential 
energy Lagrangian L = T — V , with the corresponding geodesic form 



dt 



L^=0 



or 



x 



(4) 



2 Briefly, the Euclidean SE(3)-group is denned as a semidirect (noncommutative) product (denoted by >) of 3D 
rotations and 3D translations: SE(3) := SO(3) > K 3 . Its most important subgroups are the following: 



Subgroup 



Definition 



50(3), group of rotations 
in 3D (a spherical joint) 



Set of all proper orthogonal 
3 X 3 — rotational matrices 



SE(2), special Euclidean group 
in 2D (all planar motions) 



Set of all 3 X 3 — matrices: 
cos 9 sin 8 r x 
— sin 8 cos 8 r y 
1 



SO(2), group of rotations in 2D 
subgroup of SE(2)— group 
(a revolute joint) 



Set of all proper orthogonal 
2 X 2 — rotational matrices 
included in SE(2) — group 



, group of translations in 3D 
(all spatial displacements) 



Euclidean 3D vector space 
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where subscripts denote partial derivatives, while Tj k are the Christoffel symbols of the affine 
Levi-Civita connection of the biomechanical manifold M. 

The corresponding momentum phase-space P = T*M provides a natural symplectic struc- 
ture that can be defined as follows. As the biomechanical configuration space M is a smooth 
n— manifold, we can pick local coordinates {dx 1 , dx n } £ M. Then {dx 1 , dx n } defines a basis 
of the cotangent space T*M, and by writing 9 € T*M as 9 = Pidx*, we get local coordinates 
{x 1 , ...,x n ,pi, on T*M. We can now define the canonical symplectic form wonP = T*M 

as: 

lu = dpi A dx 1 , 

where 'A' denotes the wedge or exterior product of exterior differential formsH This 2— form u> is 
obviously independent of the choice of coordinates {x 1 , ...,x n } and independent of the base point 
{x 1 , ...,x n ,p\, —,p n } G T*M. Therefore, it is locally constant, and so dco = 00 

If (P, Lu) is a 2nD symplectic manifold then about each point x £ P there are local coordinates 
{x 1 , ...,x n ,pi, ...,p n } such that ui — dpiAdx 1 . These coordinates are called canonical or symplectic. 
By the Darboux theorem, lo is constant in this local chart, i.e., dw = 0. 

3 Autonomous Hamiltonian Biomechanics 

We develop autonomous Hamiltonian biomechanics on the configuration biomechanical manifold 
M in three steps, following the standard symplectic geometry prescription (see [U [31 [H [7] ) : 

Step A Find a symplectic momentum phase-space (P,uj). 

3 R,ecall that an exterior differential form a of order p (or, a p—form a) on a base manifold X is a section of the 

p 

exterior product bundle AT X — > X. It has the following expression in local coordinates on X 

a = ctx 1 ...\ p dx Xl A ■ ■ ■ A dx x P (such that \a\=p), 

where summation is performed over all ordered collections (Ai, A p ). Cl p (X) is the vector space of p— forms on a 
biomechanical manifold X. In particular, the 1— forms are called the Pfaffian forms. 

4 The canonical 1— form 9 on T* M is the unique 1— form with the property that, for any 1— form f} which is a 
section of T*M we have (3*9 = 8. 

Let / : M — » M be a diffcomorphism. Then T* f preserves the canonical 1— form 9 on T* M , i.e., (T* f)*9 = 9. 
Thus T* f is symplectic diffeomorphism. 
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Recall that a symplectic structure on a smooth manifold M is a nondegenerate closed^ 2— form 
u) on M, i.e., for each x £ M, u)(x) is nondegenerate, and duj = 0. 

Let T*M be a cotangent space to M at m. The cotangent bundle T*M represents a union 
(J me mT*M, together with the standard topology on T*M and a natural smooth manifold structure, 
the dimension of which is twice the dimension of M. A 1— form 8 on M represents a section 
6:M->T*M of the cotangent bundle T*M. 

P = T*M is our momentum phase-space. On P there is a nondegenerate symplectic 2— form 
w is defined in local joint coordinates x l ,pi G U, U open in P, as w = dx l A dpi. In that case the 
coordinates x l ,pi € £/ are called canonical. In a usual procedure the canonical 1— form 9 is first 
defined as — p{dx % , and then the canonical 2-form lu is defined as uj = —d6. 

A symplectic phase-space manifold is a pair (P, w). 

Step B Find a Hamiltonian vector-field Xh on (P, w). 

Let (P, be a symplectic manifold. A vector-field X : P — > TP is called Hamiltonian if 
there is a smooth function F : P — ► R such that ix u) — dF (ixv denotes the interior product or 
contraction of the vector-field X and the 2-form ui). X is locally Hamiltonian if ixu) is closed. 

Let the smooth real-valued Hamiltonian function H : P — > R, representing the total biome- 
chanical energy H(x,p) — T(p) + V(x) (T and V denote kinetic and potential energy of the 
system, respectively), be given in local canonical coordinates x l ,Pi € U, U open in P. The Hamil- 
tonian vector-field Xjj, condition by ixn 1 - = dH, is actually defined via symplectic matrix J, in 
a local chart U, as 

X H = JVH = (d Pl H, -8 xi H) , J = Q , (5) 

where 7 denotes the n x n identity matrix and V is the gradient operator. 

Step C Find a Hamiltonian phase-flow <p t of Xh- 

Let (P, u>) be a symplectic phase-space manifold and Xh = JVH a Hamiltonian vector-field 
corresponding to a smooth real- valued Hamiltonian function H : P — > R, on it. If a unique one- 
parameter group of diffeomorphisms <f> t : P — > P exists so that J^|t = o0 t x = J\7H(x), it is called 
the Hamiltonian phase-flow. 

5 A p— form /3 on a smooth manifold M is called closed if its exterior derivative d = didx 1 is equal to zero, 

dp = 0. 

From this condition one can see that the closed form (the kernel of the exterior derivative operator d) is con- 
served quantity. Therefore, closed p— forms possess certain invariant properties, physically corresponding to the 
conservation laws. 

Also, a p— form that is an exterior derivative of some (p — 1)— form a, 

P = da, 

is called exact (the image of the exterior derivative operator d). By Poincare lemma, exact forms prove to be closed 
automatically, 

dp = d(da) = 0. 

Since d? = 0, every exact form is closed. The converse is only partially true, by Poincare lemma: every closed 
form is locally exact. 

Technically, this means that given a closed p— form a £ il p (U), denned on an open set U of a smooth manifold 
M any point m £ U has a neighborhood on which there exists a (p — 1)— form p S (l p ~ (U) such that dp = oc\xj. 
In particular, there is a Poincare lemma for contractible manifolds: Any closed form on a smoothly contractible 
manifold is exact. 
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A smooth curve t i— » [x t {t) 1 Pi{t)) on (P, ui) represents an integral curve of the Hamiltonian 
vector-field Xh — JVP, if in the local canonical coordinates x l ,pi £ U, U open in P, Hamiltonian 
canonical equations hold: 

q l = d Pz II, Pi = -d q iH. (6) 

An integral curve is said to be maximal if it is not a restriction of an integral curve defined on 
a larger interval of R. It follows from the standard theorem on the existence and uniqueness of the 
solution of a system of ODEs with smooth r.h.s, that if the manifold (P, uj) is Hausdorff, then for 
any point x — (x l ,pi) £ U, U open in P, there exists a maximal integral curve of Xh = JVP, 
passing for t = 0, through point x. In case Xh is complete, i.e., Xh is C p and (P, u>) is compact, 
the maximal integral curve of Xh is the Hamiltonian phase-flow <fr t : U —> U. 

The phase-flow <f) t is symplectic if w is constant along <j) tl i.e., 0^ u> = to 

{4>* t oj denotes the pull-bacf^ of lo by <f> t ), 

i&£x H oJ = 

(£x.h w denotes the L?e derivative^ of a; upon X# ). 

Symplectic phase-flow <fi t consists of canonical transformations on (P, to), i.e., diffeomorphisms 
in canonical coordinates x' l ,Pi E P, {/ open on all (P,u>) which leave w invariant. In this case the 
Liouville theorem is valid: <j> t preserves the phase volume on (P, w). Also, the system's total energy 
H is conserved along <ft t , i.e., H o <f> t = <fi t . 

Recall that the Riemannian metrics g — <, > on the configuration manifold M is a positive- 
definite quadratic form g : TM — > K, in local coordinates a;* E P, P open in M, given by (|3|) above. 
Given the metrics gij, the system's Hamiltonian function represents a momentum p-dependcnt 
quadratic form H : T*M — > M - the system's kinetic energy P(p) = T(p) = | < p,p >, in local 
canonical coordinates x l ,pi E P p , U p open in T*M, given by 

H(p) = ^g i3 {x,m)p lPj , (7) 

where g l ^(x,m) = gT^ (x,m) denotes the inverse (contravariant) material metric tensor 

Hi \ r ^ dxJ 

g J (x,m) = > m y o rs — — — — . 

X=l 

6 Given a map / : X — > X' between the two manifolds, the pullback on X of a form a on X' by / is denoted by 
f*ct. The pullback satisfies the relations 

/*(« A #) = /*« A /'ft dfa = f*(da), 

for any two forms a, (3 a f2 p (Jf). 

7 The Lie derivative SUiCt of p— form o along a vector-field « is defined by Cartan's 'magic' formula (sec 3 4 ): 

SluCt = u\dct + d(u\ct). 

It satisfies the Leibnitz relation 

£ u {a A P) = £ u a A + a A £„/3. 

Here, the contraction J of a vector-field u = u^d^ and a p— form o = a\ l ...x p dx 1 A — A cirpp on a biomechanical 
manifold X is given in local coordinates on X by 

u\a = « t 'a M A 1 ...A p _ 1 ^ Al A ■ ■ • A dx x p- 1 . 

It satisfies the following relation 

uJ(qA/3) = uJaA,3 + (-l) |a| aA«J/3. 
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T*M is an orientable manifold, admitting the standard volume form 



, N(N+1) 

n - ^ 2 ^ 

Wff ~ N~\ 

For Hamiltonian vector-field, X# on M, there is a base integral curve 7 (i) = (x z (t), Pi(t)) iff 
7o(i) is a geodesic, given by the one-form force equation 

p l =p l + Y) k gi l g km PWm = Q, with x k =g H Pi . (8) 

The l.h.s pi of the covariant momentum equation (jSJ represents the intrinsic or Bianchi covariant 
derivative of the momentum with respect to time t. Basic relation pi — defines the parallel 
transport on T N , the simplest form of human-motion dynamics. In that case Hamiltonian vector- 
field Xh is called the geodesic spray and its phase-flow is called the geodesic flow. 

For Earthly dynamics in the gravitational potential field V : M — > R, the Hamiltonian H : 
T*M extends into potential form 

H{p,x) = ^g tJ PiPj + V(x), 

with Hamiltonian vector-field Xh — JVH still defined by canonical equations ©. 
A general form of a driven, non-conservative Hamiltonian equations reads: 

x l = d Pi H, pi=Fi-d xi H, (9) 

where Fi — Fi(t, x,p) represent any kind of joint-driving covariant torques, including active neuro- 
muscular-like controls, as functions of time, angles and momenta, as well as passive dissipative and 
elastic joint torques. In the covariant momentum formulation l[8"]). the non-conservative Hamilto- 
nian equations © become 

Pi = ^ + T) k g 3l g km pip m = Fi, with x k = g kl Pi . 

The general form of autonomous Hamiltonian biomechanics is given by dissipative, driven 
Hamiltonian equations on T*M: 

dH dR 

x = ~5 ^ o - J ( 10 ) 

dpi dpi 

^ dH dR 

^(0) = 4, Pi(0)=p°i, (12) 

including contravariant equation (|10|) - the velocity vector-field, and covariant equation pip - 
the force 1-form (field), together with initial joint angles and momenta (|12[) . Here R = R{x,p) 
denotes the Raileigh nonlinear (biquadratic) dissipation function, and Fi — Fi(t, x,p) are covariant 
driving torques of equivalent muscular actuators, resembling muscular excitation and contraction 
dynamics in rotational form. The velocity vector-field (JTUJ) and the force 1— form Ijlip together 
define the generalized Hamiltonian vector-field Xh; the Hamiltonian energy function H — H(x,p) 
is its generating function. 



S 



As a Lie group, the biomechanical configuration manifold M — • SE(3y is Hausdorffd There- 
fore, for x = (x % , p{) £ E/p, where I7 P is an open coordinate chart in T*M, there exists a unique 
one-parameter group of diffeomorphisms tf> f : T*M — > T*M, that is the autonomous Hamiltonian 
phase-flow: 

<P t : T*M-.T*M:(p(0),a;(0))H-.( P (t), a; (t)), (13) 
(0t ° 4 = 0t+«) = identity), 
given by (|10H12p such that 

^ | t=0 <j> t x = JWH{x). 

4 Time— Dependent Hamiltonian Biomechanics 

In this section we develop time-dependent Hamiltonian biomechanics. For this, we first need to 
extend our autonomous Hamiltonian machinery, using the general concepts of bundles, jets and 
connections. 

4.1 Biomechanical Bundle 

While standard autonomous Lagrangian biomechanics is developed on the configuration manifold 
X, the time-dependent biomechanics necessarily includes also the real time axis K, so we have 
an extended configuration manifold Ixl. Slightly more generally, the fundamental geometrical 
structure is the so-called configuration bundle n : X — > M. Time-dependent biomechanics is thus 
formally developed either on the extended configuration manifold M x X, or on the configuration 
bundle tt : X — > R, using the concept of jets, which are based on the idea of higher-order tangency, 
or higher-order contact, at some designated point (i.e., certain anatomical joint) on a biomechanical 
configuration manifold X. 

In general, tangent and cotangent bundles, TM and T*M, of a smooth manifold M, are special 
cases of a more general geometrical object called fibre bundle, denoted tt : Y — > X, where the word 
fiber V of a map it : Y — ■> X is the preimage tt~ 1 (x) of an element x G X. It is a space which 
locally looks like a product of two spaces (similarly as a manifold locally looks like Euclidean space), 
but may possess a different global structure. To get a visual intuition behind this fundamental 
geometrical concept, we can say that a fibre bundle Y is a homeomorphic generalization of a 
product space X x V (see Figure [lj, where X and V are called the base and the fibre, respectively. 
tt : Y — > X is called the projection, Y x — n^ 1 (x) denotes a fibre over a point x of the base X, while 
the map / = tt^ 1 : X Y defines the cross-section, producing the graph (x, f(x)) in the bundle 
Y (e.g., in case of a tangent bundle, f = x represents a velocity vector-field). 

More generally, a biomechanical configuration bundle, tt : Y — > X, is a locally trivial fibred 
(or, projection) manifold over the base X. It is endowed with an atlas of fibred bundle coordinates 
{x^,y l ), where (x x ) are coordinates of X. 

A linear connection f on a biomechanical bundle Y — > X is given in local coordinates on Y by 

[H 

f = dx x ®[d x -r i jX (x)y j d i \. (14) 

8 That is, for every pair of points xi,X2 £ M, there are disjoint open subsets (charts) Ui,U2 C M such that 
xi 6 U\ and X2 £ f/2- 
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bundle Y*XxV 




Figure 1: A sketch of a locally trivial fibre bundle Y « X x V as a generalization of a product 
space I x V; left - main components; right - a few details (see text for explanation). 

An affine connection T on a biomechanical bundle Y — > X is given in local coordinates on Y 

by 

r = dx x ® p A + (-r, A (xy + T\(x))di}. 

Clearly, a linear connection T is a special case of an affine connection T. 

Every connection T on a biomechanical bundle Y — > X defines a system of first-order differ- 
ential equations on Y , given by the local coordinate relations 

y\ = r l (y). (15) 

Integral sections for T are local solutions of (fT5|) . 

4.2 Biomechanical Jets 

A pair of smooth manifold maps, /i,/2 : M — ► N (see Figure [2]), are said to be k— tangent (or 
tangent of order k, or have a fcth order contact) at a point x on a domain manifold M, denoted 
by h - / 2 , iff 

fi(x) = fi{x) called — tangent, 
d x fi(x) = d x f 2 (x), called 1 - tangent, 
d xx fi(x) = d xx f 2 {x), called 2 - tangent, 
etc. to the order k 

In this way defined k—tangency is an equivalence relation. 

A k—jet (or, a jet of order k), denoted by j£f, of a smooth map / : M — > ./V at a point ieM 
(see Figure [2]), is defined as an equivalence class of fc— tangent maps at a;, 

j'^/ : M -> N = {/':/' is k - tangent to / at x}. 

For example, consider a simple function f : X —> Y, x ^ y = f{x), mapping the A— axis into 
the Y — axis in R 2 . At a chosen point x G X we have: 
a 0— jet is a graph: (a;, /(x)); 
a 1— jet is a triple: (x, /(x), /'(x)); 
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At 

Figure 2: An intuitive geometrical picture behind the A; —jet concept, based on the idea of a higher- 
order tangency (or, higher-order contact). 

a 2— jet is a quadruple: (x, f(x), f'(x), f"(x)), 

and so on, up to the order k (where f'(x) = d ^}^ , etc). 

The set of all k— jets from j£f : X — » Y is called the k— jet manifold J k (X,Y). 

Formally, given a biomechanical bundle Y — ► X, its first-order jet manifold J X Y comprises the 
set of equivalence classes j^s, x S X, of sections s : X — ► Y so that sections s and s' belong to the 
same class iff 

Ts \t x x= Ts' | Tx x ■ 

Intuitively, sections s,s' € j^s are identified by their values s'(x) = s' l (x) and the values of their 
partial derivatives d^s l (x) — d fl s n (x) at the point x of X. There are the natural fibrations [14] 

7T i : J X Y 3 jls h-> x e X, 7T i : J X Y 3 j^s h-> s(x) € Y. 

Given bundle coordinates (x x ,y l ) of Y, the associated jet manifold J X Y is endowed with the 
adapted coordinates 

{x\y\yi), (y\y l x)Uls) = (s l (x),d,s\x)), y '\ = _ x (9 M + ^ i y i . 

In particular, given the biomechanical configuration bundle M —* R over the time axis R, 
the 1—jet space J 1 (R, M) is the set of equivalence classes j}s of sections s 1 : K — > M of the 
configuration bundle Af — > R, which are identified by their values s z (t), as well as by the values 
of their partial derivatives dts 1 = dts' l (t) at time points t G R. The 1—jet manifold J 1 (R, M) is 
coordinated by (t } x l } x l ), that is by (time, coordinates and velocities) at every active human joint, 
so the 1-jets are local joint coordinate maps 

j\s : R — ► M, t i-y (t,x l ,i l ). 

The repeated jet manifold J 1 J 1 Y is defined to be the jet manifold of the bundle J X Y — ► X. It 
is endowed with the adapted coordinates (x x 1 y l , y\, y\ p yy\^)- 

The second-order jet manifold J 2 Y of a bundle Y — > X is the subbundle of J Y — ► J X Y 
defined by the coordinate conditions y l Xfl — y^ A . It has the local coordinates (x x ,y l ,y\,y l x< ^) 
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together with the transition functions [2] 

dx a 

The second-order jet manifold J 2 Y of Y comprises the equivalence classes j^s of sections s of 
Y — > X such that 

y\(f x s) = dxs l (x), yl^s) = d^(x). 

In other words, two sections s, s' € j^s are identified by their values and the values of their first 
and second-order derivatives at the point x E X. 

In particular, given the biomechanical configuration bundle M — > R over the time axis R, 
the 2— jet space J 2 (R, M) is the set of equivalence classes jfs of sections s 1 : R — + M of the 
configuration bundle w : M — » R, which are identified by their values s l (t), as well as the values 
of their first and second partial derivatives, d t s l = d t s l {t) and d u s l = d u s l (t), respectively, at 
time points t € R. The 2-jet manifold J 2 (R, M) is coordinated by (t ), that is by (time, 

coordinates, velocities and accelerations) at every active human joint, so the 2-jets are local joint 
coordinate map^ 

fis:R^M, 1 1-» (t,x\x\x l ). 



4.3 Polysymplectic Dynamics 

Let Y — ► X be a biomechanical bundle with local coordinates y l ). In jet terms, a first-order 
Lagrangian is defined (through the standard Lagrangian density C) to be a horizontal density 
L = Cuj on the jet manifold J X Y. The jet manifold J X Y plays the role of the finite-dimensional 
configuration space of sections of the bundle Y — > X. 

Lagrangian L yields the Legendre morphism L : J X Y — ► LI, from the 1— jet manifold J l Y to 
the Legendre manifold, given by a product |14j 

n — 1 n 

n = V*Y A( A T*X) = V*Y A{AT*X)®TX, (16) 

where V is called the vertical bundle, n plays the role of the finite-dimensional phase-space of 
sections of Y — > X. Given the bundle coordinates (x x ,y 1 ) on Y — ► X, the Legendre bundle 
(116[) has local coordinates (x x ,y' l ,p x ), where p\ are the holonomic coordinates with the transition 
functions 

Relative to these coordinates, the Legendre morphism L reads 

P?oL = *V. (18) 
The Legendre manifold II in (fTH)) is equipped with the generalized Liouville form 

e = -p^dy l Auj<g>d x , (19) 

9 For more technical details on jet spaces with their physical applications, see 1 141 115| ). 
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where <S> denotes the standard tensor product. Furthermore, for any Pfaffian form 9 on X we have 
the relation 

A\e = -d(e\o). 

This relation introduces the the canonical polysymplectic form on the Legendre manifold n, 

A = dp$ /\dy i Aw® d x , (20) 

whose coordinate expression (|20[) is maintained under holonomic coordinate transformations (|17|) . 
It is a pullback- valued form [13] . 

4.4 Hamiltonian Connections 

Let J 1 !! be the jet manifold of the Legendre bundle II — > X. It is endowed with the adapted 
coordinates (x x , y % ,p x , y^Vip)- By analogy with the notion of an autonomous Hamiltonian vector- 
field Xh in (O, a connection 7 on the bundle II — > X, given by 

7 = ^0(^+7^+7^), 

is said to be locally Hamiltonian if the exterior form 7J A is closed and Hamiltonian if the form 
7JA is exact. A connection 7 is locally Hamiltonian iff it obeys the conditions [T31 [T5] 

91^-^ = 0, 9^-^ = 0, 9i7l+9l7^=0. (21) 

A p— form iJ on the Legendre bundle n — > X is called a general Hamiltonian form if there 
exists a Hamiltonian connection such that 

7JA = dH. 

A general, dissipative, time-dependent Hamiltonian form H on n is said to be Hamiltonian if it 
has the splitting 

H = H T -H T = p^dy 1 A ui\ — {p*T\ + H r )u; = p^dy 1 A lu\ — Tiui (22) 

modulo closed forms, where F is a connection on Y and Hy is a horizontal density. This splitting 
is preserved under the holonomic coordinate transformations (|17p . 

Let a Hamiltonian connection 7 associated with a Hamiltonian form H have an integral section 
s of n — * X, that is, 70 s = J x s. Then s satisfies the system of first-order Hamiltonian equations 

on n, 

Wl = W j& = -d{H. 

4.5 Time— Dependent Dissipative Hamiltonian Dynamics 

We can now formulate the time-dependent biomechanics as an n = 1— reduction of polysymplectic 
dynamics described above. The biomechanical phase space is the Legendre manifold n, endowed 
with the holonomic coordinates (t,y l ,pi) with the transition functions 

/ dy 3 
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II admits the canonical form A (|2T))) . which now reads 

A = dpi A dy l f\dt® d t . 

As a particular case of the polysymplectic machinery of the previous section, we say that a con- 
nection 

1 = dt®{d t + 1 l d l + ll d l ) 

on the bundle II — ► X is locally Hamiltonian if the exterior form 7J A is closed and Hamiltonian if 
the form 7J A is exact. A connection 7 is locally Hamiltonian iff it obeys the conditions (j2"Tj) which 
now take the form 

ay _ ay = o, d aj - d 3ll = o, erf + ay = o. 

Note that every connection F = dt <g> (dt + F l di) on the bundle Y — > X gives rise to the 
Hamiltonian connection F on n — ► X, given by 

r = dt®(dt + r i di-d j r i Pi d j ). 

The corresponding Hamiltonian form Hy is given by 

Hy = Pidy 1 - piVdt. 

Let H be a dissipative Hamiltonian form (|22p on n, which reads: 

H = Pl dy l - Hdt = p t dy l - Pl T l dt - H r dt. (23) 

We call TL and TL in the decomposition (|23p the Hamiltonian and the Hamiltonian function re- 
spectively. Let 7 be a Hamiltonian connection on n — > X associated with the Hamiltonian form 
([2"5]l. It satisfies the relations [II [T5] 

7 J A = dpt A dy' + j. t dy l A dt - j l dpi hdt = dH, 

y = a ,; ^, 7j = (24) 

From equations (|24[) we see that, in the case of biomechanics, one and only one Hamiltonian 
connection is associated with a given Hamiltonian form. 

Every connection 7 on n — > X yields the system of first-order differential equations (fl~5|) , which 
now takes the explicit form: 

f=l\ Pi = li- (25) 
They are called the evolution equations. If 7 is a Hamiltonian connection associated with the 
Hamiltonian form H (1231) . the evolution equations (|25p become the dissipative time- dependent 
Hamiltonian equations: 

f = d l H, Pi = -diH. (26) 

In addition, given any scalar function / on n, we have the dissipative Hamiltonian evolution 
equation 

d H f = (d t + d l Hd % - dm 1 ) f, (27) 
relative to the Hamiltonian H. On solutions s of the Hamiltonian equations (f2"6"l) , the evolution 
equation IpT)) is equal to the total time derivative of the function /: 

s*d H f=^-(fos). 
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4.6 Time & Fitness— Dependent Biomechanics 

The dissipative Hamiltonian system ([26]) - ([27f is the basis for our time & fitness-dependent biome- 
chanics. The scalar function / in (|2"T)) on the biomechanical Legendre phase-space manifold II is now 
interpreted as an individual neuro-muscular fitness function. This fitness function is a 'determi- 
nant' for the performance of muscular drives for the driven, dissipative Hamiltonian biomechanics. 
These muscular drives, for all active DOF, are given by time & fitness-dependent Pfaffian form: 
Fi = Fi(t,y,p, /). In this way, we obtain our final model for time & fitness-dependent Hamiltonian 
biomechanics: 

if = &H, 
pi = Fi - diH, 
d H f = (8 t + d l Hd i -dSHd i )f. 

Physiologically, the active muscular drives Fi = Fi(t,y,p, /) consist of Q][2]) : 

1. Synovial joint mechanics, giving the first stabilizing effect to the conservative skeleton 
dynamics, is described by the (y,y)-form of the Rayleigh-Van der Pol's dissipation function 

n 

i=l 

where on and (i i denote dissipation parameters. Its partial derivatives give rise to the viscous- 
damping torques and forces in the joints 

P° int = dR/dy\ 

which arc linear in y l and quadratic in y % . 

2. Muscular mechanics, giving the driving torques and forces F™ usc = Tl ausc {t,y,y) with 
(i = 1, . . . , n) for human biomechanics, describes the internal excitation and contraction dynamics 
of equivalent muscular actuators [llj . 

(a) The excitation dynamics can be described by an impulse force-time relation 
pvra-p = pQ(i_ e -t/r^ if stimulation > 
F un P = p o e -t/r, if stimulation =0, 

where F® denote the maximal isometric muscular torques and forces, while Tj denote the associated 
time characteristics of particular muscular actuators. This relation represents a solution of the 
Wilkie's muscular active-state element equation [12] 

fi + Tfi = TSA, fj,(0) = 0, 0<5<1, 

where /i = /j,(t) represents the active state of the muscle, T denotes the element gain, A corresponds 
to the maximum tension the element can develop, and S = S(r) is the 'desired' active state as a 
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function of the motor unit stimulus rate r. This is the basis for biomechanical force controller. 

(b) The contraction dynamics has classically been described by the Hill's hyperbolic force- 
velocity relation [13] 

(SijV 3 + h) ' 

where a? and bi denote the Hill's parameters, corresponding to the energy dissipated during the 
contraction and the phosphagenic energy conversion rate, respectively, while 5ij is the Kronecker's 
S— tensor. 

In this way, human biomechanics describes the excitation/ contract ion dynamics for the ith 
equivalent muscle-joint actuator, using the simple impulse-hyperbolic product relation 

JT usc {^v,y) = Fr p x F " m - 

Now, for the purpose of biomedical engineering and rehabilitation, human biomechanics has 
developed the so-called hybrid rotational actuator. It includes, along with muscular and viscous 
forces, the D.C. motor drives, as used in robotics 

Tl oho = » fc (t) - J k y k (t) - B k y k (t), with 
lkh{t) + R k ik(t) + C k y k (t) = Ufc(t), 

where k — 1, . . . , n, i k (t) and u k (t) denote currents and voltages in the rotors of the drives, R k , l k 
and C k are resistances, inductances and capacitances in the rotors, respectively, while J k and B k 
correspond to inertia moments and viscous dampings of the drives, respectively. 

Finally, to make the model more realistic, we need to add some stochastic torques and forces: 

Tf° ch = B tJ [y t (t),t]dWHt), 

where Hy [y(t) , t] represents continuous stochastic diffusion fluctuations, and W J (t) is an TV— variable 
Wiener process (i.e., generalized Brownian motion) [5], with 

dW j (t) = W°{t + dt) - W j (t), (for j = l,...,n = no. of active DOF). 



5 Conclusion 

We have presented the time-dependent generalization of an 'ordinary' autonomous Hamiltonian 
biomechanics, in which total mechanical + biochemical energy is not conserved. Starting with 
the Covariant Force Law, we have first developed autonomous Hamiltonian biomechanics. Then 
we have introduced a general framework for time-dependent Hamiltonian biomechanics in terms 
jets, Legendre manifolds and dissipative Hamiltonian connections associated with the extended 
musculo-skeletal configuration manifold, called the configuration bundle. In this way we formulated 
a general Hamiltonian model for time & fitness-dependent human biomechanics. 
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